Effect of electronic interactions on the persistent current in one-dimensional 

disordered rings 
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The persistent current is here studied in one-dimensional disordered rings that contain interacting 
electrons. We used the density matrix renormalization group algorithms in order to compute the 
stiffness, a measure that gives the magnitude of the persistent currents as a function of the boundary 
conditions for different sets of both interaction and disorder characteristics. In contrast to its 
non-interacting value, an increase in the stiffness parameter was observed for systems at and off 
half- filling for weak interactions and non-zero disorders. Within the strong interaction limit, the 
decrease in stiffness depends on the filling and an analytical approach is developed to recover the 
observed behaviors. This is required in order to understand its mechanisms. Finally, the study of 
the localization length confirms the enhancement of the persistent current for moderate interactions 
when disorders are present at half-filling. Our results reveal two different regimes, one for weak and 
one for strong interactions at and off half-filling. 

PACS numbers: 73.23-b,73.23.Ra,71.10.Fd 



I. INTRODUCTION 

Mesoscopic rings are known to carry an equilib- 
rium current called persistent currents when they are 
threaded through a magnetic flusi. The experimental 
discoverjiSi^ of persistent currents has motivated a num- 
ber of theoretical studies. For example, non-interacting 
electron theories^^Sii describe the general features but 
they strongly underestimate the current present within 
the diffusive regime*. The fact of taking into account ci- 
ther electron-electron interactions'^ or simply the degree 
of disordei"^ leads to an underestimation and hence, only 
a partial description of the magnitude of the currents. A 
possible way to improve the characterization of the per- 
sistent currents is to take into account both disorder and 
electronic interactions-^'^. Unfortunately, this has proved 
to be a difficult problem to tackle. Diagrammatic pertur- 
bative expansion in three dimension suggests an increase 
in the current.'^ but does not provide the means to con- 
clude firmly. Therefore, we restricted the problem here 
to one dimension. Surprisingly, the role of electronic in- 
teractions is still a controversial question and some have 
even suggested that such interactions may lead to a de- 
crease in the magnitude of the persistent currentsi2*i^. 
On the one hand, experiments using exact diagonaliza- 
tion for systems containing very few lattice sitesi^ and 
the self-consistent Hartree-Fock approximationMii& show 
that repulsive interactions suppress the persistent cur- 
rents whereas calculations show no dramatic increase in 
the magnitude of the persistent currentsi^iii. On the 
other hand, DMRG numerical studies carried out on spin- 
less fermions claim to obtain a small increase in the mag- 
nitude of the persistent currents for very high disorder 
only^*. It has also been shown that repulsive interactions 
may provoke significant disorder and hence, enhance the 
electron mobilitjiiSiSS in few-particle models. The need 
to take into account the spin has also theoretically been 



investigated using the method of bosonization~. The re- 
sult suggests that the stiffness could be enhanced by in- 
termediate interactions. Numerical approaches^SiS^ fur- 
ther show that repulsive interactions enhance the persis- 
tent currents in weakly disordered systems off half-filling 
only. 

In this paper, we consider a simplified model of in- 
teracting electrons in a one-dimensional disordered ring, 
which is penetrated by a magnetic flux. It has been re- 
ported in another study that moderate interactions could 
induce an enhancement of the persistent currents even 
in the presence of weak disorder for half-filled systems 
Moreover disorder could have the unexpected effect of in- 
creasing the zero-temperature persistent currents in one- 
dimensional half-filled Hubbard rings, with strongly in- 
teracting electrons. We emphasize that even off half- 
filling the persistent current may still be enhanced by 
small interactions and in such case, the filling may have 
an influence on its behavior. Indeed, we show that the 
behavior of the persistent currents changes as it is moved 
off half-filling, owing to the appearance of holes. These 
results are established partially within the framework of 
the strong coupling perturbation expansion and partially 
through numerical calculations by means of the Density 
Matrix Renormalization Group algorithni2SiS&. Finally, 
we study the influence of the spin on the behavior of the 
stiffness parameter and show that the case of half-filling 
is rather particular. This manuscript^^ is organized as 
follows. In Section im the model and the calculation al- 
gorithms of the persistent currents are presented. The 
numerical and analytical work at half-filling are presented 
in section Hill Section ITVI contains the numerical results 
off half-filling and further explains the results by means 
of the perturbation expansion. Section Ivl highlights the 
influence of the spin not only on the stiffness parame- 
ter, but also on the magnitude of the persistent currents. 
Section IVII is devoted to the study of the localization 
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length and identifies two regimes that are characterized 
by weak and strong interactions, which are separated by 
a crossover phase. Finally, our conclusions are presented 
in section rvnl Details of the analytical approaches used 
in Sections IIIII and IIVI are given in Appendices lA II and 
lA 21 A toy model is developed for non half-filled systems 
in Appendix lA 31 Appendix ^ concerns the behavior of 
the weights of the ground-state in the strong interaction 
limit. In Appendix[ni the specific case of an odd number 
of electrons present in a ring is discussed. 



II. 



MODEL AND METHOD 



A one-dimensional lattice is shown in function of the 
Hubbard- Anderson HamiltonianSS-, which we use to de- 
scribe N interacting electrons in a disordered ring of M 
sites: 
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denotes the kinetic energy, t will be set equal to one in 
the numerical part. 
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stands for the disorder with Vi e [^^, W > and 
finally 
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(4) 



accounts for the electronic interactions. 

The threading of the ring by a magnetic flux introduces 
an Aharonov-Bohm phase $ = 2n(j)/(f)() measured in units 
of the flux quantum (po = he/ e. According to Bayers 
and YangSS the only effect of the magnetic field is to 
impose twisted boundary conditions co,cr — CA/,crexpi$ 
between the first and the last site, the other pieces of the 
Hamiltonian remain unchanged. The persistent currents 
that are carried through an isolated ring are the direct 
consequence of the sensitivity of the eigenenergies to a 
phase-change in the boundary conditions, in particular 
for 



m 



(5) 



at vanishing temperature where i?o('/') denotes the 
ground-state energy. It has been shown that transport 
properties are related to thermodynamic properties in 
a localized system at vanishing temperature'^''. In fact, 
persistent current are a periodic function of the flux 0, 



/($) = (£;(0) - £;(7r))sin$ 



(6) 



where 



E{0) - E{tt) = AE 



(7) 



is called phase sensitivity. We furthermore introduce the 
stiffness, given by 



D 



M. 



\AE\ 



(8) 



which is related to the Drude weight Dc . We compute 
D because it is simpler to calculate than the conductivity 
CT, and because it furthermore allows for a measure of the 
magnitude of the persistent currents. The energy levels 
are computed by means of the density matrix renormal- 
ization group method and exact diagonaliza- 

tion. The first version of the program was written by P. 
Brune^i for the ionic Hubbard model. Since H is invari- 
ant under rotation around the z-axis its eigenstates can 
be calculated by restricting it to the subspace — Q- We 
used eight finite-lattice iterations and kept a maximum 
of 700 states per block iteration. 



III. 



HALF-FILLING 



The case of half-filling is here considered when the 
number of sites is equal to the number of electrons 
M ^N. 

The very strong interaction limit can be treated analyt- 
ically in order to better comprehend the behavior of the 
phase sensitivity within the Mott insulator limit. When 
U ^ W, t the interaction U dominates, and it is possi- 
ble to expand AE in terms of t/U^'^-^^. The details of 
the calculation are given in Appendix lA II The basis- 
states I'^a) of Hq spanning the a-space'^^ are given by 
the product of the different on-site states, and read as 
follows: 




|0), 



The function il}-^\a) specifies the site where the k 
electron with spin t (i) among the N/2 t (i)-electrons 
is created following the configuration a. The |0) state 
is composed of M empty sites. The symbol j^'/j) de- 
notes states that do not possess doubly-occupied sites 
(exactly one electron per site is present at half-filling), 
spanning the 5-space; l^'^y) denotes states that contain 
at least one doubly-occupied site spanning the I?-space. 
The corresponding many-body energies of the states {ipa) 
are summed over the disorder contributions E^^ and the 
interaction energies E^. They are equal to: 
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Here ga stands for the number of doubly-occupied sites. 
In the Mott insulator hmit when U :^ W,t the |5'-y)'s 
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possess a very high energy level and hence, a gap arises 
with the |*/3) s. 

First, we consider the case without disorder. The 
|\E'/3)'s are the basis-states and all possess an equal zero 
energy level. Hence, for very large U, H — Hq + Hw and 
the kinetic energy term Hk represents a perturbation 
with respect to the dominating term Hij. Within the 
limit oi U ^ t, the spin degree of freedom can be treated 
as a significant spin Hamiltonian that is in fact the an- 
tiferromagnetic Heisenberg Hamiltonian when (/) = 0, for 
the case of the clean systemS^iSSi^. The antiferromag- 
netic Heisenberg Hamiltonian ground-state is the super- 
position of all states with a total spin equal to zero. It is 
furthermore non-degenerate. This theorem has been ex- 
tended to the half-filled Hubbard model"^^-^^ for an even 
number of sites M^Si^ and the ground-state can be ex- 
pressed as: 




|0) (11) 
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where the < fp < 1, according to Marshall's theorem' 
for all spin configurations f3 in the subspace S and i^, ^ 
i\.,yk ^ k' . The ground-state for <i> = tt is equal to that 
presented in eq. ((TT|l . 

Now we calculate the leading contributions to the 
phase sensitivity Ai? using the difference present in the 
higher-order corrections of the many-body ground-state 
energies under periodic and antiperiodic boundary con- 
ditions. Studying the Hubbard Hamiltonian within the 
strong interaction limit, is equivalent to studying the ac- 
tion of Hk in the 5-spacc'^*. Indeed, the numerator of 
is composed of matrix elements (^^7, IH/fl "071+1) of 
the hopping part of the Hamiltonian, which connects 
each site to its neighbor. Thus, the resulting correc- 
tion brought to the ground-state energy is different 
from zero only when all 1^7) states can be connected by 
one-particle hopping processes^Si^. These sequences are 
called A^'^'^ ^ and start at and end at \'4'i3'), with 
the intermediate states being all \ilj^ys. In order to cal- 
culate the difference between periodic and antiperiodic 
boundary conditions, the particles must circle the ring 
and cross the border once. The lowest contributions are 
of order M, and in such case, the energy corrections for 
periodic and antiperiodic conditions differ only by the 
sign. Intermediate results are given in an already pub- 
lished article24. 

Second, we reintroduce disorder. The 5-states follow 
Mott-insulator configurations and all possess a similar 
and lowest energy level that reads Eq^ = W X]f=T^ '^i 
(from eq. HlOfl 'l. This energy level is independent of 
the spin configurations of the 5-states and consequently, 
within this limit, they degenerate and are characterized 
by the same weights fi}{U 00). All the terms in the 
denominator are calculated by differences between dis- 
order and interaction contributions of the ground-state 
l^'o) and the intermediate state |4'7, ). The denomi- 
nators are written as Wd^^ + 9-yiU. We further define 



{E^ — E^) — Wd~^^, as the energy difference between 
the disorder configurations of the state |\l/o) and the state 



We develop AE in powers of jj- up to second or- 



der. This yields : 
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This sum is performed over all forward hopping sequences 
^(/3,/3 ) ^•£gj. electrons that are moving in the anti- 
trigonometric direction of motion) with a first factor 
2 introduced for the backward sequences and an addi- 
tional factor 2 accounting for the difference between pe- 
riodic and antiperiodic contributions. All contributing 
sequences lead to a cyclic perturbation of the N/2 opera- 
tors corresponding to electrons with a specific spin direc- 
tion. Since the weights fp are all positive, the sign of the 
phase sensitivity at strong interaction is given by (— 1)^^^ 
in eq. (|12|l . as in the non-interacting case^. The domi- 
nating term in eq. (|12|l (corresponding to the clean limit 
W 0) exhibits an interaction-induced suppression of 
the persistent current cx U{t/U)^'^. The first-order cor- 
rection in W/U totally vanishes due to the particle- hole 
symmetry. The second-order correction is positive be- 
cause the first term is positive, and much larger than the 
second one. We conclude that the phase sensitivity has 
positive corrections for order (^)^. This conclusion is 
also true for the stiffness parameter, as D = ^ | Ai5| . The 
unexpected conclusion is that at half-filling, the disorder 
increases both the phase sensitivity and the magnitude of 
the persistent currents for highly interacting electrons in 
Hubbard- Anders on rings^'^. This counterintuitive effect^ 
might also be related to a competition between the phe- 
nomenon, disorder and interaction. The ground-state en- 
ergy, which is equal to £^0 = ''^i increases with the 
strength of the disorder W, and when W < U, the energy 
level becomes comparable to the interaction strength U. 
Under these conditions, the gap reduces between the en- 
ergy levels of the excited states j^^) and the ground- 
state. A jump to a doubly-occupied site is then rendered 
easier and thus, the movement of the electrons around 
the ring is favored, enhancing the magnitude of the per- 
sistent currents. These sets of M hops illustrate a simple 
mechanism that can explain the changes in the magni- 
tude of the persistent currents under strong interaction 
limits. 

Extensive numerical calculations were performed for 
systems of iV = 14 electrons (7 t; 7 i) on M = 14 sites, 
and for a sample of = 20 particles (10 t; 10 i) on 
M = 20 sites. 

A selection of the probability distributions of D{U) are 
plotted for fixed M and N values, for set disorder and 
interaction values. They can be fitted by Gaussian func- 
tions and are said to be log-normal. In Fig. ^ the vari- 
ance cr^ is plotted as a function of U, for M = N = 14. 
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FIG. 1; Variance of the logarithm of the stiffness as a function 
of U, for W = 2,5,9 and for M = N = 14. 



For U — 0, the fluctuations are governed by the disorder 
and they increase when disorder increases. When inter- 
actions are introduced, fluctuations of two different na- 
tures compete. Figure ^ illustrates this by showing that 
these fluctuations compensate each other and decrease, 
phenomenon that leads to an increase in the magnitude 
of the persistent currents. Overall, this produces a max- 
imum of magnitude for the persistent currents around 
U ^ W/2 ± t. For moderate U, the particles can move 
and two electrons can be together on the same site. Thus, 
the electrons explore sites with higher on-site potentials 
and under such conditions, the fluctuations are more sen- 
sitive to disorder. Finally, in the region where U < W, 
the variance increases. After this second maximum, the 
fluctuations decrease when the interaction increases for 
U > W. Within the region of very strong interactions, 
the particles are trapped alone within their sites and thus, 
their movement is frozen with exactly one particle per 
site. Fluctuations arc on-site only and decrease system- 
atically. 

Fig. Elshows the dependency of the stiffness parameter 
log 13 on the degree of interaction, for M = = 14,20 
for typical individual samples. These samples are repre- 
sented by the single lines and global means are illustrated 
for three different disorder values, W, and for the clean 
case where W = (solid line). The ensemble averages 
of log D are performed over approximately 100 different 
samples (disorder realizations) . The behavior of both the 
samples and the means are similar but in contradiction 
with the spinless fermiongS. Indeed, our curves are all 
rather smoothSii^ even if small peaks appear for the in- 
dividual samples for strong disorders {W ^ 7, 9) in Fig. 
121 This second peak increase of \ogD for samples oc- 
curs for U ^ W — t when the fluctuations are strong 
(cf. Fig. For U = 0, the disorder leads to the An- 
derson localization and logD is strongly suppressed by 
increasing disorder—. This is what would be expected in 
one dimension, for a non-interacting disordered system. 
For clean rings {W = 0), the interaction always reduces 
the stiffness, which is consistent with a Luttinger liquid 
calculationMi. 



FIG. 2; Dependency of the stiffness parameter upon the de- 
gree of interaction, for M = iV = 14 (M = iV = 20 in the 
inset), for disorders W = 0,2, 5, 9. Lines represent individual 
samples. Symbols correspond to the global average. Error 
bars are smaller than the symbols size and hence, do not ap- 
pear. Note that the curve for W = behaves as a Luttinger 
liquid. 



A weak repulsive interaction U > leads to an in- 
crease of the stiffness compared to its value for J7 = 0, 
even when little disorder is present; such an increase was 
in fact predicted from a renormalization group approach, 
for both off half-filling and moderate disordersSi. Indeed, 
the stiffness is less sensitive to disorder when repulsive in- 
teractions are accounted for than when attractive inter- 
actions only are considered, off half filling and for mod- 
erate disorders. This is clearly illustrated in Fig. O for 
an individual sample M = A^ = 10. For attractive inter- 
actions no increase occurs and the stiffness phenomenon 
is strongly suppressed. The curves corresponding to dif- 
ferent disorders do not cross to negative interactions, in 
accordance with the analytical development of Sec. IIIII 
which is valid only for positive interactions"^^. At inter- 
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FIG. 3: Dependency of the stiffness parameter for an individ- 
ual sample upon the degree of interaction U , for M = A'^ = 10 
and for disorders W = 0,2, 5, 9. Interactions are both attrac- 
tive and repulsive, —20 < U < 25. An increase occurs for 
positive interactions only. 
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mediate interaction levels i.e. for U G VF], the phase 
sensitivity exhibits a maximum around W/2 ± which 
becomes broader with increasing disorder. When dis- 
order increases this maximum decreases, but the ratio 
Draax/D{Q) is greater. With M = iV = 14, the mean 
ratios are equal to 1.4, 8.1 and 150 for W = 2,5 and 9, 
respectively. These ratios are even greater for samples 
with M — N = 20 and can under certain conditions be 
three times greater (e.g. for = 9). 

For large repulsive interaction U 3> W , the behav- 
ior of \ogD is radically different. The stiffness decreases 
strongly with the degree of interaction and the numer- 
ical results confirm this with the analytical power law 
(X l/t/*^~^, as illustrated in Fig. 0| The slopes for mod- 
erate interactions are stronger when the disorder strength 
W increases in order to recover the l/C/*^~^ behavior 
when J7 ^ VF. Results for clean systems {W = Q) and 



M=N=li ■ 



M=N=20 




1.1 logioU 



FIG. 4: Data are plotted on log-log scale for M = = 14 
{M = Ai' = 20 in the inset) under strong interaction condi- 
tions. Results confirm the power law D <x f/-(A^-i). Symbols 
illustrate global means and lines individual samples. 



strong interactions gave the means to check a good con- 
cordance between observed results and the power law. 
This decrease is related to the localization due to the 
interaction^. The most surprising result is that a curve 
crossing occurs within the Mott insulator limit (see Fig. 
121) confirming that the disorder induces an increase in the 
stiffness parameter, as found analytically. 

For all values of U and W, we obtained negative AE for 
M ^ N ^ 6, 10, 14, and positive AE for M = N = 20. 
Thus, consistent with both Leggett's rule^ for the non- 
interacting case and the analytical result for very strong 
interactions, the sign of AE is determined by the number 
of electrons that are each characterized by a given and 
specific spin, according to (—1)^/^. Our numerical re- 
sults revealed that this rule is valid for all tested samples 
and for all considered interaction and disorder values. In 
the following section, these results are compared to the 
phenomenon observed off half-filling. 



IV. OFF HALF-FILLING 

This section concerns systems off half-filling, N ^ M . 
The very strong interaction limit can be treated analyt- 
ically and we maintain the notations presented in Sec. 
mil Details of the calculations are given in Appendix 
lA 21 An example is developed in a toy model in Ap- 
pendix El The expansion is valid for an even number 
of particles and a single empty site, called a hole^SiSS. 

Once again, we first consider the case without disorder. 
For very large U the kinetic energy term Hk represents 
a perturbation with respect to the dominating term Hjj 
as at half-filling^. We must first determine the ground- 
state of the system. The basis-states are still given by eq. 
@. Now since a hole has been introduced, the 5-states 
are connected by only one hop. A given basis-state \^ p) 
can be connected to itself within an ensemble of first or- 
der hops, which are all going in the same direction to form 
a super-lattice"*^. For the toy model of two electrons on 
three sites of Appendix I A 31 the number of super-lattices 
is equal to one. Increasing both the number of sites and 
the number of electrons increases the number of super- 
lattices and leads to the creation of N/2 subblocks, which 
is the number of electrons with the same spin present 
within the ring. We call these subblocks Ci, where i is 
the number of electrons of similar spin ones besides the 
others. Each subblock d leads to a different eigenvector, 
which implies that the ground-state degenerates. The 
basic structure of the ground-state can be determined by 
the standard Perron-Frobenius sign conventionifi. in the 
case of a single empty site N — M — \. The ground- 
state is in such case a superposition of the 1 5'^' ) within 
each subblock Ct . The sign of the weights depends on the 
boundary conditions^SiS. If the transition amplitude be- 
tween two states is negative, then one can superimpose 
within a subblock, states with the similar signs. Vicev- 
ersa, if the transition is positive, one can superimpose 
states with opposite signs. These are in fact the Marshall 
sign rules^ that are extended to the Hubbard model—. 
The ground-state inside each subblock is written as: 




> for $ = 
e for $ 



IT 



(13) 



where the weights fp are chosen to be real^L^. Thus, 
for antiperiodic boundary conditions within each sub- 
block d, half of the basis-states will possess a positive 
(or negative) sign^SiSS. The first order corrections of the 
many-body ground-state energies are equal for periodic 
and antiperiodic boundary conditions. Indeed, in case of 
antiperiodic boundary conditions, each time an electron 
crosses the border, the sign due to the flux is compen- 
sated for by the sign of the weights. The phase sensitivity 
cancels out at first order. At second order, the intermedi- 
ate states are obviously equal to |^'7)'s, with exactly one 
doubly-occupied site. Lets now consider those sequences 
that cross the border. For example, that sequence where 
a particle jumps to form a doubly-occupied site and then. 
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one of these electrons conies back to the initial site or 
jumps onto a neighboring site. These sequences are called 
exchange and double-jump. Their relative terms, when 
considering the two Hamiltonians, have opposite signs 
for periodic and antiperiodic boundary conditions. The 
eigenvalues for these matrices are different and hence, the 
phase sensitivity reads: 



A£;(2) = |^(2)(0)-S(2)(^) 



(14) 



where the factor ^j- appears in all parts of both matrices. 
The phase sensitivity behaves as 1/U within the very 
strong interactions limit. 

Now, if one introduces disorder, then the particles stay 
on the sites with the lowest disorder potentials inside each 
subblock Ci and form the C/*-subspace. But one main- 
tains N/2 subblocks. Hence, the ground-state is now a 
superposition of the basis-states with the weights fp 
respecting the Marshall rulc'^^'^^. A single jump does not 
lead to a basis-state anymore and one considers it now 
as two jumps. Under such conditions, the only possible 
sequences are the exchanges that lead to opposite values 
for periodic and antiperiodic boundary conditions. The 
phase sensitivity at second order here reads, after an ex- 
pansion of W/U : 



(15) 



This 1/U power law is lower than at half- filling. Indeed, 
for N — M, the ground-state is rigid. When N < M, 
with the provided hole, the electrons can move more eas- 
ily around the ring, which decreases the power law com- 
pared to the half-filled situation. This nicely illustrates 
the behavior of the persistent currents off half-filling un- 
der strong interaction limits. Electrons move around the 
ring without meeting each other and when close to the 
border, they hop on a doubly-occupied site. The phase 
sensitivity then decreases as the disorder strength W in- 
creases. Unfortunately, the perturbation theory does not 
allow us here to conclude about the sign of AE under 
strong interaction limits. 

We now present numerical simulations for = 8 par- 
ticles (4 t; 4 I) on M = 9 and on M = 11 sites. The sym- 
bols stand for the mean performance calculated over 100 
disorder realizations (within a given sample); the lines 
represent individual samples. The results concentrate on 
an even number of electrons N, and we specifically con- 
sider the effect of the parity of the number of sites M. 
For an odd and an even value of M , the behavior of the 
stiffness is similar. Moreover, samples present similar be- 
haviors as that observed for the mean curves. Finally, all 
curves are relatively smooth. The averages are calculated 
over the logarithm of the stiffness, \ogD. 

The variances for three values of the disorder W — 
2,5,9 are plotted as a function of the interaction for 
iV = 8 particles, for M = 9 and M = 11 sites. Re- 
sults show that they are similar in Fig. |S1 For < ?7 < 
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FIG. 5: Variance of the logarithm of the stiffness log D as a 
function of U for different strengths of the disorder W = 2,5,9 
for M = 11 (M = 9 in the inset) and N = 8. 



W/2 ± t, the behavior is similar to that expose in Sec. 
mil with the fluctuations that decrease. The persistent 
currents reach a maximum around W/2 ± t. Considering 
the strong interaction limits, the particles are trapped 
alone within their sites but a movement is still permitted 
because of the existence of an empty site (hole). The sys- 
tem is thus not rigid. Fluctuations are more significant 
than at half-filling since the fluctuations of both nature 
do not compensate, i.e. the on-site fluctuations and the 
fluctuations that are due to the disorder present on the 
empty site. When the disorder increases the fluctuations 
in the empty site increase, and so does a^. 

In Fig. El the logarithm of the stiffness log D is plot- 
ted as a function of the interaction U, for different val- 
ues of the disorder W and for the clean case specifically. 
At U — \ogD is strongly suppressed by the disor- 




FIG. 6: Dependency of the stiffness upon the degree of in- 
teraction U, for N = 8, M = 11 (M = 9 in the inset) and 
W = 0,2,5,9. Lines represent log D for individual samples 
and symbols correspond to global means. 

der (Anderson localization)^. For clean rings, we verify 
that the interaction systematically reduces the stiffness, 
which is consistent with a Luttinger liquid calculation^. 
A weak repulsive interaction leads to an increase of the 
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stiffness respectively to its value in absence of interac- 
tions, as predicted by a renormalization group for mod- 
erate disordersSi. We obtain this increase for all fillings 
and all non-zero disorder strengths. This increase can 
reach factors equal to approximately 1.5 (1.3), 4.5 (3.4), 
28 (17) and when W = 2,5 and 9, respectively. This is 
true for the average system M = 11 (M = 9) and N = 8. 
The increases of the samples are important and can reach 
an order of magnitude. 

The stiffness log I? decreases as the interaction U in- 
creases, within the very strong interaction limits. Figure 
[3 presents the logarithm of the stiffness log Z? as a func- 
tion of log U for M = 11, iV = 8 and gives the means 
to determine that the power law follows 1/U. We have 
verified that for the system M = 9,N — 8, the stiffness 
decreases as a function of C/~^, in agreement with strong 
interaction situations. For important disorder conditions, 
the slope must be calculated for interactions U ^ W in 
order to recover the l/U law. This law is true for sev- 
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FIG. 7: Data are plotted for strong interaction conditions in 
log-log scale for M = 11 and N = 8. They reveal a power law 
of 1/17. The symbols represent means and the lines, individual 
samples. 

eral holes even if the perturbation theory is valid only 
for single holes. In Fig. |S1 the stiffness is smaller for 
negative interaction values than for positive onesSi, and 
no stiffness increase occurs. 

The sign of the phase sensitivity was also studied and 
results demonstrate that it also follows the Leggett's 
rules^. The sign depends on the number of electrons 
that possess a specific spin, according to (— 1)^^^. This 
was observed for all disorder and interaction values. The 
present result was verified numerically for all the com- 
puted samples. For N = 8, the phase sensitivity is always 
positive and for = 6 the phase sensitivity is negative, 
whatever the number M of sites. 



V. ROLE OF THE SPIN 

The logarithm of D gave us the means to model the 
behavior of the persistent currents as a function of U, 
for different values of the disorder W and for different 
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FIG. 8: Illustration of the behavior of log D for a sample of 
M = 9 sites and N = 8 electrons for negative and positive 
interactions —20 < U < 20. Note that an increase occurs for 
positive interactions only. 



system sizes M. Our conclusion is that a pronounced 
increase occurs for all degrees of disorder not equal to 
zero, for every system size M and arbitrary fillings, but 
for moderate values of the interaction U ~ W/2. 

This increase can be explained through the competi- 
tion between disorder and interaction phenomenon. In- 
deed, the disorder sets the electrons within the sites 
with the lowest potentials whereas the interaction sets 
a unique electron per site. When disorder and interac- 
tion are of the same strength, the determination of the 
ground-state of the system is difficult. Indeed, the elec- 
trons oscillate between being together on the same site 
hence, adding an energy level to U or being on two differ- 
ent sites which sets a different degree of disorder, energy 
Wvi. However, the ground-states directly depend on the 
spin parameter, where the case of half-filling corresponds 
to M = A^. When W > U, electrons may prefer avoiding 
the site with highest potentials and hence, will generate a 
doubly-occupied site (a spin up and a spin down). When 
U and W are of the same order, the electrons move eas- 
ily as the doubly-occupied site can hop freely around the 
ring and set itself next to the site with the highest poten- 
tial. Under such conditions, ii W ^ U, a. particle from 
the doubly-occupied site can jump on the site and fur- 
ther on to its neighbor, recreating a doubly-occupied site. 
The quantity here important to consider is the difference, 
U — W/2. When it is smaller than the amplitude tran- 
sition t, the electrons movement around the ring is easy 
and induces an important current. This behavior can 
be related to calculations for 2D spinless fermions with 
disorder and long range Coulomb interaction^^. When 
W < U, the particles remain alone on their respective 
sites and this leads to a decrease in the magnitude of the 
current. For off half-filling, this explanation is also valid 
but the holes are placed upon the highest disorder po- 
tentials. Hence, the maximum current further decreases. 
This simple mechanism of two competing ground-states 
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is a starting point for a better understanding of what 
happens in two dimensions, and provides a starting in- 
dication of how persistent currents may be enhanced by 
large orders of magnitudes. 

We illustrate this modeling considering first a very 
weak disorder W = 0.1, for the two following systems 
M = TV = 10, and = 8, M = 11. Here logZ? is cal- 
culated from a group average of about 100 samples. An 
increase is already observable for such a weak disorder in 
Fig. 121 One can see that a maximum increase is obtained 
for U ~ 0.05 = W/2. These results give us the means to 



0.5 



0.4 



0.3 



♦ M=11,N=8 
M=N=10 


^ 

♦ 

♦ ♦♦♦ ^ * 






<> 





♦ 













U 



FIG. 9: Interaction dependency of the stiffness parameter on 
the degree of interaction, for W — 0.1 in the case of two 
systems, M = iV = 10 and iV = 8,M = 11. These results are 
calculated from a group average of 100 samples. 

conclude to an increase of the stiffness for non- vanishing 
disorders. The persistent currents are less important for 
M = 11, N = 8 than for M = N = 10. 

For negative interactions, the stiffness always de- 
creases. Indeed, two electrons with opposite spin direc- 
tions can both generate a doubly-occupied site and be 
on the sites with lowest potentials. This decreases the 
stiffness as soon as [/ < for all disorders strengths, as 
found in Fig. O 

In Fig. ^1 a sample of iV = 4 electrons on M = 13 
sites is presented in order to show that even for low 
fillings i.e. [N/M] ^ 0.3, persistent currents increase. 
This increase reaches an order of magnitude for moder- 
ate {W = 5) and strong {W = 9) degrees of disorder. 
Thus, the mechanism that have been presented above 
to explain the increase of the persistent current is valid 
even for low fillings and for attractive and repulsive in- 
teractions in presence of disorder. Other studies have 
previously considered the role of the spin and have also 
concluded that it is an important parameter for both one 
and two dimension situations, and for multi-channelled 
systemsi^. The present work shows in addition that in 
one dimension, the spin has a significant influence on the 
magnitude of the persistent currents. 

In the previous sections Ijllll and llV|l . the limit for the 
strong interactions revealed different behaviors for half- 
filled and non half-filled systems. We now consider the 




FIG. 10: log D is plotted as a function of U , for different 
values of the degree of disorder W , for an individual sample 
of = 4 electrons present on M = 13 sites. 
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FIG. 11: logD is illustrated as a function of U, for W = 
{W = 5 in the inset) for systems with N = 8 and M = 9, 10 
or 11. A significant crossing between the curves is observed 
for strong interactions. 



case of a similar number of particles and a similar de- 
gree of disorder (even W = 0) but with different filling 
properties (i.e. not equal to 0.5). For strong interactions 
i.e. when the filling decreases for a given number of elec- 
trons, the stiffness increases. This is illustrated in Fig. 
111! The presence of holes makes the persistent currents 
increase within the strong interaction limit. The case of 
half-filling deserves to be treated apart and presents a 
different behavior, which underlines the important effect 
of the particle-hole symmetry for the elcctrona^^. 



VI. 



LOCALIZATION LENGTH 



In this final section, we study the localization length ^ 
in order to characterize the finite-size effecli^. In order 
to extract ^, we use a relation that is valid for localized 
systems. The localization length is related to the stiffness 
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D through: 



\nD = \n{A{U)) - M/S,{U). 



(16) 



\ogD is plotted as a function of M and one must check 
that a straight hne is obtained. Thanks to a hnear re- 
gression the slope, which is proportional to can 
be obtained from the numerical data of log!?, indicated 
here by the dimension of the error bars. This relation is 
phenomenological and is obtained thanks to dimensional 
reasoning but does not allow for the determination of the 
sign. 



A. Localization length at half-Ailing 

Disordered non-interacting one-dimensional systems 
are Anderson insulators^^, whereas a half- filled sys- 
tem behaves as a Mott insulator for very strong 
interactions^. We consider half-filled rings M — N = 
6, 10, 14. We first verified the results obtained for W — Q. 
First, plots representing \ogD{M) must be straight. For 
[/ = = 0, a free electron gas is obtained and the local- 
ization length should both diverge and be infinite. Be- 
cause of the numerical errors, one gets a very high finite 
value of the localization length ^. Once the interaction 
is introduced, ^ can be extracted. It decreases as a func- 
tion of the interaction increase (see Fig. I12|) . This has 
previously been reported^^. For moderate interactions. 




FIG. 12: (, is plotted as a function of U, for W = Q. The 
symbols represent the values of ^ obtained from logD. In the 
bottom inset, f A— is presented as a function of the interac- 
tion U , for 3 < f7 < 8. In the top inset, ,f ~^ behaves according 
to the Bethe Ansatz within the strong interaction limit^. 



8*2 

U 



the 



(7 e [3,8], the quantity eA (A- 

Lieb-Wu charge gap^^) must reach the mean-field result 
of {£,A)/t — 4*^. The values calculated thanks to our nu- 
merical approach reached 4.9, for U — (see the bottom 
inset of Fig. For very strong interactions the in- 

verse of the localization length behaves as the analytical 
solution given by the Bethe Ansat»4i. This is in agree- 
ment with what was presented in Sec. IIIII The inverse 



of the localization length ^ ^ (see the top inset of Fig. 
[Hjl behaves as -1.49 + lnU for W = 0. 

We want to study the influence of the competition be- 
tween disorder and interaction on the localization length 
^. The inverse of ^ is represented in Fig. E|as a func- 
tion of \nU, for different strengths of the disorder W. 
The inset presents log i:»(M) for U = 1,W = 2. The 
U = limit shows that ^ decreases when the disorder 
increases as for the Anderson localization. For moderate 
interactions, an enhancement of the localization length ^ 
when compared to its non-interacting value is observed 
for all disorder strengths W. The values for W = 2 are 
very high, even larger than the sizes of the studied sys- 
tems. In such a case, eq. ()16f) is not valid anymore. 
Yet, it is possible to conclude to a strong and signifi- 
cant increase of ^. The increase factors ^inax/^(C^ — 0) 
reaches 2.2, 3.4 and 4.5, for W = 2,5 and 9, respec- 
tively. They are significant and confirm the delocaliz- 
ing effect. The maximum of £,{U,W) belongs to the 
interval ± t] and decreases as W increases. This 
strong enhancement suggests that the wave functions are 
less localized and that the electrons might move more 
easily on several sites. Moreover, the stiffness depends 
on the localization length according to the equation: 
D{M) = {~l)^^/^+^M^/^D{U,W)exp-M/£_{U,W)M.. 
Hence, when ^ increases, D also increases and the per- 
sistent currents will consequently increase even for large 
finite rings of circumference M. Another weak increase 
for the localization length had been reported for spinless 
fermions but only for a strong degree of disorder—. In 
the case of strong interactions, the localization length ^ 
decreases with increasing J7, as a Mott insulator. As for 
the stiffness D that was presented in Sec. IIIII close to 
the strong interaction limits, the disorder phenomenon 
has the unexpected effect of increasing the localization 
length. This phenomena that is observed for systems of 
AI — 20 sites can be extended to bigger systems. 

In order to illustrate both competing regimes, one 
should fit the inverse of the localization length. Fig. 
1131 shows that depends on \nU/t for all disorder 
strength. This dependence is expected for very strong 
interactions and W = 0^. However, for W and for 
weak interactions, the inverse of the localization length 
behaves like lnJ7~l°l as a correlation length whereas for 
strong interactions it behaves following InU^^^^. For ex- 
ample when W — 2, behaves as —1.0 -I- O.SlnC/ 
for weak interactions and then, as 0.022 — 0.001 \nU for 
strong interactions, in Fig. ^| In Fig. ^1 log 13 is rep- 
resented as a function of In U and two different regimes 
also clearly appear. First D increases up to C/+l° ' and 
then decreases up to [/+''' L This is true for all the sys- 
tems that we considered as well as for all non-zero dis- 
orders. The energy density e — dE/dM was represented 
as a function of the interaction U , for different disor- 
der strengths W = 2,5,9 (see inset of Fig. ^J. At 
half-filling, when U is strong, the particles tend to stay 
isolated within their site and the system becomes rigid, 
as previously mentioned in Sec. IIIII The energies con- 
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interaction reduces, electrons can jump on other sites and 
the crossover can then occur. For weak interactions, the 
electrons can move easily on multiple sites and the move- 
ment is dominated by quantum fluctuations. The second 
derivative of the quantities logi? and e illustrate 

this crossover. For weak interactions, d^e/dU^ behaves 
as — C^exp— C?7 while for strong U, it is a constant. 
Between these limits, a peak exists and symbolizes the 
crossover. A study of the shift for this peak as a func- 
tion of the degree of disorder should be considered to 
describe the corresponding crossover as a function of the 
interaction, Ucmss as shown in Fig. 1181 



FIG. 13: C^{M = N) is plotted in function of \nU . The 
dotted lines represent the fitting of both regimes for W — 2,5 
and 9. 
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FIG. 14: logD (e) is plotted as a function of InU (U) for 
half- filled systems for W = 2, 5, 9. The dotted lines are the 
fitting for both regimes. 



verge towards zero, and the density of energy increases 
as — exp— CC/ ~ CU. For very strong interactions i.e. 
when there is exactly one particle per site, there is lit- 
tle displacement and the energies correspond to the total 
sum of the on-site potentials. As the systems are half- 
filled and the potentials are random, the sum tends to 
converge to zero for all degrees of disorder. The energy 
density tends to saturate at zero. Thus, one fits e for 
weak and strong interactions with straight lines propor- 
tional to U. 

Figures |13I and mi present the two competing regimes 
and reveal a very sharp crossover, as the fittings joining at 
Ucmss covers almost the totally of the symbols used in the 
different curves. Here, we illustrate two exclusion limits. 
For U — 0, the electrons with opposite spin remain on the 
same site (Pauli principle) while for infinite interactions, 
the electrons avoid being on the same site (the exclu- 
sion principle). For strong interactions, the mobility of 
the electrons is weak since they are trapped alone within 
their site and hence, the movement is frozen. When the 



B. Localization length off half- filling 

The study of the localization length is now considered 
for off half-filling. One considers three systems with a 
similar number of particles N = 8, but with different 
size M = 9, 10, 11. Some of these systems were presented 
in Sec. |IV| 

Fig. 1151 illustrates log D{M) for different degrees of 
interaction, with W = 5. As in Fig. ^] results reveal 
that for U > W, log D{M) increases as M increases. The 
slope is thus positive. This leads to a change in the sign 
of the localization length off half-filling which is in agree- 
ment with the comments made following eq. (|16l) . One 




FIG. 15: logD is plotted as a function of M — 9, 10, 11 for a 
degree of disorder W = 5, for N = 8 particles. The slope of 
log D changes sign as U increases. 

can assume that approximately straight lines are present 
even if only three points are here considered. At the point 
where the curves logD{U, M) cross (see Fig. Illll . the 
localization length strongly increases and reaches a high 
positive value. After the crossing, it reaches an important 
negative value before increasing once more. Figure [TCI il- 
lustrates this behavior for W = 0,2, 5, 9, considering the 
inverse of the localization length When U = W = 0, 
the localization length starts at a very high value and de- 
creases when U increases. When U = 0, ^ decreases as 
the degree of the disorder increases as it is case for the 



11 



Anderson insulator—. For W = 0, ^ decreases, then, a 
change in sign occurs and ^ increases up to a disorder- 
dependent negative value. For all disorders W ^ and 
moderate values of U, ^ increases strongly compared to 
its non-interacting value. If one considers the positive 
values of ^, it is noticeable that ^(0) > C(2) > ^(5) > ^(9) 
(except for [/ = 2, due to numerical errors). If the high- 
est positive value is not considered, the ratios 6nax/C(0) 
reach 1.5, 2.75, and 4 for = 2,5 and 9, respectively. 
These ratios are close to but smaller than those observed 
at half-filling. After the change in sign, the localization 
length increases up to a negative value ^nog that depends 
on the degree of disorder. When the disorder increases, 
the localization length decreases as expected for off half- 
filhng (Sec. ITVIl . 

In order to study both regimes (Pauli and exclusion 
limits) as was in Sec. IVI Al one fits as a function of 
In U for weak and strong interactions in order to obtain 
the J7cross- This is presented in Fig. ^1 For positive 
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FIG. 16: is plotted as a function of in U ior W ^ 0, 2, 5, 9, 
for N = 8 particles. The straight lines correspond to the fits 
of both regimes. A crossover is here identified. 

values, the inverse of the localization length behaves 
as InJ/"!"^! and for strong interactions, the localization 
length behaves as lnJ7~l''^l because of the change in sign. 
In Fig. El the curves plotting log 13 as a function of 
In U, for various degrees of disorders, are represented 
and it is possible to fit them with lines proportional to 
In U. They present behavior similar to that observed for 
half-filling behaviors. The energy density e — dE/dN 
is represented as a function of the interaction U , in the 
inset of Fig. El For weak interactions, the electrons are 
placed by pairs (with opposite spin directions) within the 
same site. In addition, the corresponding energies are 
summed over iV/2 sites. When the interactions become 
stronger, the energies increase as it would be the case for 
sums over less negative potentials. The energies of the 
ground-state when the number of sites increases, decrease 
strongly. This leads to a decrease in the energy density 
as expressed by exp— CJ7 ~ —CU. For very strong in- 
teractions, the ground-state energies are summed over 
on-site potentials. As there are empty sites, the energy 
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FIG. 17: log D (e) is plotted as a function of In U (U) for 
N — 8. The dotted lines correspond to the fitting of both 
studied regimes. 

density tends to saturate at a specific negative value that 
increases with the strength of the disorder. 

The quantities e, D and f highlight the presence of a 
crossover as mentioned in Sec. IVI Al and gives the means 
to determine the existence of a peak between the two dif- 
ferent regimes. The position of the peak is given by the 
value C/cross, when the curves fitting both regimes cross. 
In Fig. El the movement of the peaks is a function of 
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FIG. 18: Ucross as a function of the disorder W for systems 
at and off half-filling. 

the strength of the disorder VF, for systems at and off 
half-filling. When W increases, the interaction must be 
stronger in order to reach the second regime where the 
particles will be localized alone within their site. The 
peaks and hence, the crossovers occur quickly around the 
same value of the interaction J7, for a given disorder W . 
This is true for all systems. The values of C/cross are very 
close but more scattered for off half-filling. The competi- 
tion between the two regimes that are represented by the 
Anderson and the Mott insulators is very important in 
order to gain a better understanding of the two dimen- 
sional metal-insulator transition, which is at the heart 
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of the problems in mesoscopic physics^^. With disorder 
only, no metal-insulator transition is foundiS.. When in- 
cluding these interactions, it is possible to determine two 
different regimes in one dimension. This may explain the 
metal-insulator transition in two dimension situation. 



VII. CONCLUSION 

We investigated the effect of electron-electron on-site 
interactions on the persistent currents found in disor- 
dered one-dimcnsional rings with arbitrary fillings. We 
took into account the spin and considered the stiffness 
as a measure of the magnitude of the persistent cur- 
rents. Our numerical simulations indicate a strong in- 
crease in the stiffness with respect to the non-interacting 
case observed for moderate interactions and arbitrary fill- 
ings. This strong enhancement, which is also observed at 
the level of disordered realization, can be attributed to 
the mechanism of both competing forces, i.e. disorder 
and interaction. A simple computation comparing the 
weakly interacting state systems in the presence of dis- 
order shows an increase in the electrons mobility within 
the ring. The stiffness decreases for stronger interactions 
according to a power law that is different at and off half- 
filling. The analytical approaches helped us to determine 
the power laws for both cases. An unexpected behavior 
was also observed at the strong interaction limit, specif- 
ically at half-filling: the disorder rendered the stiffness 
stronger. Such an effect can be explained by the disorder- 
induced reduction in the energy gap between the Mott 
insulator ground-state and the excited states. By study- 
ing the localization length, we here further observed that 
the interactions increase the persistent currents even for 
a large finite size of rings at half-filling. The localiza- 
tion length, which is defined phenomenologically, changes 
sign off half-filling. The behavior of the energy density, 
the stiffness as well as the localization length suggest a 
crossover between weak and strong interactions. Based 
on these results, we can conclude that the spin plays a 
significant role and influences the increase in magnitude 
of persistent currents. 



APPENDIX A: PERTURBATION EXPANSIONS 
1. Half-Ailing 

We give here the details for the analytical approach in 
Sec. mil The ground-state is given by eq. ((TT|) for both 
periodic and antiperiodic condition o'^^i'^^ . The aim here 
was to define sequences (the set of hops defined in Sec. 
IIII|I that may lead to discrepancies within the corrections 
brought to the energy level. 

The terms presented in the numerator of i?^") are equal 
to —t but the signs of the permutations of the elec- 
trons, with spin up or down, could be negative as the 



hopping element changed the electrons order. For peri- 
odic boundary conditions, the numerator can be writ- 
ten as Np = (— i)"'sign(P|)sign(F|) where we define 
P|(I)(A) td) as the permutation of the positions for up 
(or down) electrons around the ring, which is a resultant 
from the defined sequence A. Since the flux appears 
only within the kinetic part Hk of the Hamiltonian, a 
unique change only appears when a particle crosses the 
border for antiperiodic boundary conditions. This in- 
troduces a sign (—1) in the numerator. Thus, the nu- 
merator for antiperiodic boundary conditions is equal to 
: Nap — (—1)'*'' Np where hi, is the number of hops 
across the border 1 < — > M contained in A. For or- 
ders lower than M, the n'^ order energy is a constant in 
<i> that does not depend on the boundary conditions and 
the phase sensitivity is zero. Indeed, an even number of 
hops is needed to recover one of the eigenstates and in 
this case, the border is crossed twice. To obtain a dif- 
ference, the particles must circulate around the ring and 
the lowest contributions are of order M; in such case, a 
particle crosses the boundary condition only once. Our 
work gave us the means to define different sequences con- 
necting two basis-states, |\l//3) and |^^), and to obtain eq. 
(|12|l . These sequences have previously been published—. 

When disorder is reintroduced, the perturbation the- 
ory remains second order, with an exchange energy that 
is now equal to rrp^ . The matrix ele- 

ments of the Hamiltonian are under such conditions even- 
functions of so are the weights fp. They read: 

fp « fj^^ + fj^\^)'^- In the following development, we 

maintained /^"^ only and f/3 ~ fj^^ was approximated. 
In order to study the effect of disorder on the phase sen- 
sitivity for strong interactions {U ^ W), we expand the 
denominator under powers of W/U . When second order 
is reached, this yields : 

13,13' ^W.fi') -li'^T^' 

(-fi:^.f?(i:|.E^))(M, 

The dominating term of (|A1|I corresponds to the power 
law found in section lTTll We now consider the term of the 
first order in W/U and this first-order correction for W/U 
vanishes because of the particle-hole symmetry. Indeed, 
if one characterizes a given sequence of hoppings by the 
positions of the doubly-occupied and the empty sites, one 
can always construct moving the particle along the ring a 
second sequence by exchanging the places of the doubly- 
occupied and the empty sites. These two sequences will 
possess an equal number of doubly-occupied sites g^, and 
their coefficients will have opposite signs. This first 
order term vanishes when the sum over all the sequences 
is taken. Thus, the second order term of ljAl|l determines 
the disorder dependence of the persistent current for U 3> 
W . The first term is of order 1 and is always positive. 



13 



Hence, for each sequence the sum over positive quantities 
of order one only are taken and it is possible to estimate: 
< M. By summing over all the sequences, we 

obtain a positive term of order NsM. The second term 
is also of order 1 but has a random sign. Taking the 
sum over all sequences, it is evaluate to almost ^^^^ . 
Since the number of sequences is very high, one obtains 
VTVj <C Ns ■ This results in eq. H12I) . It is concluded here 
that the phase sensitivity possesses positive corrections 
of order (■^)^, so does the stiffness D. 



2. Off half-filling 

We give here the details of the calculations for the per- 
turbation expansion that was presented in section IIVI 
The ground-state is given by eq. ((T^ . 

The matrix element of second order contains inter- 
mediate states with exactly one doubly-occupied site. 



|'07)'s, that connect |V'^*) to IV"^"), as the second jump 
leads to l^'/j). We seek to develop sequences of two 
jumps that will lead to different terms accord- 
ing to the defined boundary conditions. Then, crossing 
border sequences are considered. The first sequences are 
the exchange sequences, where two electrons exchange 
places with one another, across the border, as illustrated 
in Fig. 1191 Under such conditions and for both peri- 
odic and antiperiodic boundary conditions, the factor 

F = {ip'^p'\HK\^p-y){ipj\HK\ipp') is equal to More 
specifically, under periodic boundary conditions, weights 
have equal signs, whereas under antiperiodic boundary 
conditions, weights have opposite signs because of the 
Marshall's rules^S*^. A second type of sequence, called 
the double jump is illustrated in Fig. ^1 and brings 
different terms Ti'*'^^ to the Hamiltonians. Under pe- 
riodic and antiperiodic boundary conditions, the factor 
F is equal to t'^ and — t^, respectively. In addition, un- 
der periodic boundary conditions, weights maintain equal 
signs. However, under antiperiodic boundary conditions, 
the double jump is equivalent to two hops for a given par- 
ticle i.e. equivalent to a change in sign and an exchange. 
This is also illustrated in Fig. ^| Overall, this sec- 
ond type of sequence gives an end product with an equal 



sign for the weights of both /' 



Ci 



and . 



This leads to 



a contribution of a second order phase-sensitivity. The 
factorization of —t^/U for all the different terms of both 
Hamiltonians leads to eq. (|14|l . 

When disorder is introduced, the ground-state is a su- 
perposition of the basis-states within each subspace 
with weights chosen to be real. At this point, one 
assumes that the Marshall rule is still valid in presence 
of disorde»224^. Due to signs of the weights, the only 
possible sequences that could lead to opposite values for 
periodic and antiperiodic boundary conditions are the 
exchange sequences. Now, using sums over two-hops se- 




FIG. 19: Example of a "double jump" sequence for M = 5 
and A'^ = 4. The third and fourth drawing illustrates an 
"exchange" sequence. 



quences A'-^'^ \ the second-order phase sensitivity reads: 



(A2) 



The expansion of eq. ljA2|l leads to the eq. (|15|l . 

3. Toy model 

The toy model is presented here to help explain the 
stiffness behavior that was presented in Sec. IIVI for gen- 
eral systems of N electrons on M sites. Here, we con- 
sider a simplified model of two electrons for three sites. 
Such a model was studied by H. Tasaki and provided the 
means to verify parts of our results^^. This small sys- 
tem is described by the Hubbard- Anderson Hamiltonian 
of eq. when restricted to three sites. In this case, 
one considers one empty site onlji^LS, as Sec. IIVI 

We first use numerical manipulations in order to com- 
pute the many-body ground-state energies and to deduce 
the phase sensitivity. The basis-states are given by eq. 
©. Thus, diagonalization of the matrices is necessary to 
obtain the eigenstates for all values of the parameters t, 
U and Vi. The case without disorder is now considered. 
Under periodic boundary conditions, the energy level is 
described by S(0) = -t- ^VlP^T^Ut + m^. As far 
as the antiperiodic conditions are concerned, the minor 
eigenvalue is equal to E{tt) = —2t. Now it is possible to 
calculate, using eq. {T)), the contributions to the phase 
sensitivity. This would lead to: 



AE 



E{0) - E{Tr) 



= t+\U - \y/{U'^ + Wt + 36t2) (A3) 



For very strong interactions, an expansion in t/U can be 
performed and the phase sensitivity is then given by: 



u 



(A4) 



This quantity behaves as l/[/ and decreases when U in- 
creases. For periodic boundary conditions, the ground- 
state of the system is a superposition of the states p) 



14 



and 1^'-,.). All weights are positive. Indeed, the unnor- 
malized weights for the states j^*^) are all equal to one 
whereas those for the states \^^) are equal to /*. One 
assumes Af^ = 3(/*)^ + 6. The components of the |'I'7)'s 
are written as follows: 



iW + t- + AtU + 36^2 



(A5) 



The development of these weights for very large U is 
such that (see eq. (^^): 



8i 
2UM' 



(A6) 



The weights of the |5'-y)'s decrease down to zero whereas 
the weights of the states are equal to 



(A7) 



We consider that for the ground-state and within the 
limits of very strong interactions, one can neglect these 
|5'-y)'s. Under antiperiodic boundary conditions, the cor- 
responding eigenvector is the superposition of the |^'/3)'s 
with alternated signs. 

Second, it is possible to calculate the corrections to 
bring to the energy level using the perturbation the- 
ory in order to recover the values of the phase sensi- 
tivity given by eq. (|A4p . The basis-states remain the 
iS-states. The first state has an t-electron on the first 
site and a |-electron on the second site, the last site be- 
ing empty in all cases. The other states are generated 
by one hop processes. This set of six states constitutes 
a " super-lattice"— and is a good illustration of the con- 
nectivity condition2&. The basic structure of the ground- 
state is determined following the Perron-Frobenius sign 
convention^Si^Sil^. The ground-states are in fact the su- 
perposition of the six states, denoted |1) • • • |6), with signs 
respecting the Marshall's rule^i as indicated in the fol- 
lowing: 

I*0,*=0(*=7t)) = ^ 

(|l) + |2) + |3> + H|4) + (-)|5) + (-)|6)) (A8) 

We can now calculate the periodic and antiperiodic 
many-body first-order energies levels. Under periodic 
boundary conditions, each |^/;^) is coupled to two other 
states by a single jump, and all terms are positive 
£■(0) = \{~t) X 2 X 6 = —2t. Under antiperiodic bound- 
ary conditions, the negative sign that is induced by the 
flux, is compensated by the negative sign of the weights 
each time the electron crosses the border line. One gets 
E{tt) = —2t. This result shows that the phase sensitivity 
cancels at first order, as explained in Sec. IIVI Similar 
eigenvalues are obtained through the numerical diagonal- 
izing of the matrices that are formed with the basis-states 
|1) • • • |6). The calculations must be increased to the sec- 
ond order in order to recover a contribution. Hence, 
one must consider two-hops sequences across the border. 



The first jump leads to a doubly-occupied site whereas 
the second jump leads to one of the many basis-states. 
Each state \i) {i — 1, ... ,6) is connected to two states 
1*-^). Each of these |^'^)'s can lead to 4 different \^p)'s. 
Under periodic boundary conditions, since the weights 
are equal and posses similar signs, the contribution is 
£;(2)(0) = ^[6x2x4x (-^2)] = -8t^/U. In the case 
of antiperiodic conditions, the sum over these terms can- 
cels out. Indeed, if we now consider the state |1), half 
of the terms contain a negative sign and the other half 
a positive one and hence, the contribution is zero. A 
similar comment can be given to other states, thus, the 
phase sensitivity is equal to —St^/U as was obtained in 
eq. (jA4|) . This confirms the validity of the perturbation 
theory within the limit for very strong interactions. 

In this toy model one introduces a disorder v on the 
first site, when considering the limits for very strong in- 
teractions. The 6*6 matrix can be diagonalized with the 
disorder v placed upon the first site. This yields : 



-Aperiodic — -^antiperiodic 



{E{v)+E{-v)) 
2 



^jt- iVw^ - 2tv + 9^2 - iV-y^ + 2to-h9t2. (A9) 

Even in the presence of disorder, the first-order phase 
sensitivity is equal to zero. Hence, for strong disorders, 
particles occupy the disordered states associated to the 
following weight: 



//3 



" W' 



9^2 



where (7V')^ = Epifp?- The ff^iW oo) decrease 
when W increases. This gives the means to avoid the 
sites with highest potentials in Sec. IIVI 



APPENDIX B: VERY STRONG INTERACTIONS 
LIMIT 

This part concerns the infinite interaction case. The 
ground-state components at half-filling were considered. 
A program was developed in order to calculate the 
weights of the eigenvectors when two particles can not 
remain together upon the same site (U — > oo). We com- 
puted the weights of each component for the ground-state 
as a function of the ring size M, for M G [2,14]. The 
most significant weights, /ait, are those of the states with 
alternated spins; the second most important ones, /max2, 
correspond to states where two electrons have exchanged 
their positions when compared to the states with alter- 
nated spins. We have represented the logarithm of these 
/alt and /max2 Weights as a function of the size of the sys- 
tem M (see Fig. I2U|I . It was established that the weights 
f{M) can be fitted by an exponential. This result gave 
the means to extrapolate what would be obtained in the 
case of large N: we speculate that for big systems, the al- 
ternated spin configurations will be the most significant 
contributor to the phase sensitivity parameter. More- 
over, as the ratio R saturates and as soon as M > 12, 
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one may conclude that the corrections brought to the 
phase sensitivity will essentially depend on both /ait and 
/max2 weights. Here, the particles move easily around 
the ring, especially for the states with alternated spins. 
Consequently, an important current may be induced. 



FIG. 21: log D as a function of the interaction U for W — 
(0, )2, 5,9 ioT M = N = 7 (for M = 11, = 7 in the inset). 
For intermediate levels of interactions, a curve increase is ob- 
served. Specifically, for the sample and W — 9, & minimum 
occurs at logD — —4.2. 




10 N 



FIG. 20: Logarithm of the most significant components for 
the ground-state of half-filled systems (A*' — M) (in the inset 
is presented the logarithm of their ratio R = /ait//max2) for 
infinite interactions as a function of the system size, M. 



APPENDIX C: ODD NUMBER OF ELECTRONS 

Here, we consider the rings filled with an odd number 
of electrons. The behavior is presented in Fig. |^ for 
two systems Af = 11, iV = 7 and A/ = iV = 7. For a 




1 ■:*i********i**ir-k 

M=ll, N=7 



disorder equal to zero, the phase sensitivity is equal to 
zero for the half-filled case, M = N = 7. The results 
show for the sample, an absence of Anderson localization 
for U = 0. Moreover the sign of the phase sensitivity de- 
pends on the on-site disorder potentials and is random. 
This renders the averaging process difhcult. Mean av- 
erages must be performed on AE, and the fluctuations 
are in such a case much more significant than for sys- 
tems with an even number of electrons. This is shown 
in Fig. 12] The averaged stiffness terms (logD) respect 
Anderson localization for J7 = for both systems. For 
moderate values of the interactions and for the disor- 
ders W = 2,5, 9, an increase of \ogD is observed for the 
mean calculated averages and some various studied sam- 
ples (except here for = 2, 5 at half- filling) . But in 
all cases, log D is smaller than that observed for an even 
number of electrons. Nevertheless logD can reach an or- 
der of magnitude (I?max/-D(C^ = 0) ~ 34), at half-filling 
(Af — N — 7) and for W — 9. For stronger interactions, 
the stiffness parameter decreases with U. At half-filling, 
the disorder increases the magnitude of the persistent 
current within the very strong interaction limits as what 
was observed for the case of even numbers of electrons. 
Finally, a perturbative approach is not possible as those 
theorems are applied uniquely for the case of even num- 
bers of both sites and electrons. 
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